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The charm contribution to the proton structure, F£ (x,Q 2 ), is determined using the in- 
clusive cross sections of D* ± (2010) meson production in deep- inelastic scattering. The 
cross section measurement covers the region 5 < Q 2 < 1000 GeV 2 in photon virtuality 
and 0.02 < y < 0.70 in the inelasticity of the scattering process. The D* ± meson is 
measured in transverse momentum and pseudo-rapidity down to p T > 1.5 GeV and up 
to \rj\ < 1.5. The data were taken with the HI detector corresponding to an integrated 
luminosity of 347 pb _1 . F% is determined from the D* production cross sections and 
compared to leading and next-to-leading order perturbative QCD predictions. 
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1) ' 1 Introduction 

The charm quark production in clcctron-proton (ep) scattering is dominated by the boson- 
C"| \ gluon-fusion (BGF) process (p/p — ► cc). The charm contribution, F^x^Q 2 ), to the proton 

structure is obtained by using the expression for the one photon exchange cross section for 
charm production: 

" (l + (l-J/) 2 ) F^Q 2 )- (1) 



> 



OO dxdQ 2 

yjQ , The y range of the measurements is such that the contribution from the longitudinal struc- 

Cn| ' ture function F£ is negligible. The current analysis [T] uses data taken with the HI detector 

qq . [2] during the HERA-II running period at a centre-of-mass energy of yfs = 318 GeV. Pre- 

liminary results pHH] of the inclusive D* cross section measurements at medium and high Q 2 
0\ . are used. Compared to early results [5] by HI a larger phase space and a significantly larger 

integrated luminosity of 347 pb _1 has been used. Hence, more precise tests of perturbative 
QCD (pQCD) become possible. 
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$_i ■ 2 Theoretical Models of Open Charm Production 

a | 

The description of heavy flavour production in ep collisions is based on pQCD at leading 
order (LO) or at next-to- leading order (NLO) for which calculations in several schemes 
are available [6l [8] - In this analysis, the D*^ meson cross sections in the visible range, as 
well as the charm contribution to the proton structure function are calculated using two 
different models. On the one hand the HVQDIS program [6] is used. It is based on a NLO 
calculation in the fixed-flavor-number-scheme (FFNS) providing differential cross sections 
of massive charm quarks including mass effects at the production threshold. The parton 
evolution is performed according to the DGLAP evolution equations [?]. HVQDIS applies 
independent fragmentation. On the other hand the LO Monte Carlo program CASCADE 
[T2], supplemented with parton showers, based on the CCFM evolution scheme is used. The 
CCFM evolution equations [13] are expected to be more appropriate to describe the parton 
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evolution at small x. The hadronisation of partons is performed using the Lund String model 
as implemented in PYTHIA [T7]. In case of HVQDIS the proton parton density function 
(PDF) MRST04FF3nlo [18] is used, while for CASCADE the proton PDF parametrisation 
AO [19] is used. The renormalisation and factorisation scale for both pQCD calculations has 
been set to fx 2 . ^ = Am 2 + Q 2 . 



2.1 Fragmentation Model 

Differential inclusive cross sections of charmed mesons are calculated using the HQVDIS 
program after fragmenting the charm quarks into D*^ mesons. The Kartvelishvili frag- 
mentation function [20j is used, which is controlled by a single parameter a. This is also 
used for CASCADE, where higher charm resonances as determined by ALEPH [5T] are in- 
cluded. The charm fragmentation function was measured at HI using inclusive D*^ meson 
production associated with jet production, where a different behaviour of the fragmentation 
function close to the threshold of charm production and far above the threshold is reported 
[22] . This results in different values for a (see Table [1), which suggests a s-dependent charm 
fragmentation function, s denotes the invariant mass of the produced cc pair. Since the 
fragmentation model influences the kinematic distributions of the D* meson an additional 
uncertainty from the fragmentation model is assigned. 



3 D** Cross Section Data & Extraction of F*(x, Q 2 ) 

The cross section measurement covers the kinematic region of 5 < Q 2 < 1000 GcV 2 and 
0.02 < y < 0.70 in the inelasticity of the scattering process. The visible range of the D*^ 
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Figure 1: The double differential cross section in y and Q 2 compared to the prediction 
from HVQDIS a) and to the prediction provided by CASCADE b). The s-dependent charm 
fragmentation function is used for both predictions. The shaded band in a) reflects the 
theoretical uncertainty estimated by parameter variations given in Table [TJ 
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meson measurement is restricted to pt{D*) > 1.5 GeV and |r/(D*)| < 1.5. The double 
differential cross section as a function of Q 2 and y is shown in Figure |TJi) compared to the 
prediction from HVQDIS and in b) for CASCADE. Both describe the data reasonably well. 
A total experimental error of 9% was achieved. These data are used for the extraction of F$. 
The measured inclusive D*^ cross sections cr™ p (y, Q 2 ) in bins of y and Q 2 are converted to 
a bin centre corrected F^ exp ((x), (Q 2 )) by the relation: 

OXp «*>, <Q 2 » = °Z^ Q nl ■ F 2 the0 «*>> (Q 2 )) > ( 2 ) 

CT vi S [y>Q ) 

where cr*j s eo and F-f tllco are the theoretical predictions from the model under consideration. 
The Bj0rken variable x is related to y via x = Q 2 /y ■ s. The contribution of open beauty 
production to the visible D*^ meson 
cross sections is estimated to be of 
the order of 1-2% and is therefore ne- 
glected. The extraction of F^ exp is 
faced with an intrinsic problem as the 
measurement covers about 30% of the 
total phase space for charm produc- 
tion. The ratio aj^/a^ is the ex- 
trapolation factor to the full phase 
space which depends strongly on the 
underlying model. The ratio of the 
extrapolation factors from CASCADE 
normalised to HVQDIS is shown in 
Figure [21 where in general an agree- 
ment between the models at the level 
of 10% is observed. This is not the 
case for the highest x region, where 
differences in extrapolation of 80% are 
observed, which is known to depend 
strongly on the r) region of the D* Figure 2: The extrapolation factors from CAS- 
cross section measurement. If ex- CADE are normalised to the ones obtained by 
tended to \<r)(D*)\ < 1.9 this differ- HVQDIS. 
ences go down to 20% [23] . 

3.1 Extrapolation Uncertainties 

The extrapolation needed to determine F% ° xp depends significantly on the underlying model 
and introduce extrapolation uncertainties. These are estimated by varying the model param- 
eters like charm mass, renormalisation and factorisation scales, PDF set and fragmentation 
parameters. The parameter variations together with the average relative uncertainties are 
summarised in Table [T] 

3.2 Results 

F£ determined from the inclusive D*^ cross sections as a function of x for different values 
of Q 2 is shown in Fig. [3^) or b) using either HVQDIS or CASCADE for the extrapolation 
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name 



variation (HVQDIS) 



rel. unc. variation (CASCADE) rel. una 



charm mass 

scale jUo = Q 2 + ^ m c 

PDF 



1.3 <m c < 1.6 GeV 
0.5 < Hrj'/no < 2 
CTEQ5f3 & MRST04ff3 



6% 
4% 
13% 



1.3 <m c < 1.6 GeV 
0.5 < fir,f/lM> < 2 
AO-, AO , A0+ 



7% 
2% 
2% 



I fragmentation: 



low s(< 70 GeV 2 ) 
high s{> 70 GeV 2 ) 




a = 3.3 ± 0.4 
70 ± 20 GeV 2 



a = 8.2 ± 1.1 
a = 4.6 ±0.6 
70 ± 20 GeV 2 



s threshold 



5 % 
2 % 



8% 
2% 



Tabic 1: Uncertainty of quantities which is utilised for the estimation of the extrapolation 
uncertainty calculated from HVQDIS and CASCADE. The relative uncertainties on the 
extrapolation factor are summarised. 

to the full phase space. The inner error bar corresponds to the statistical error, whereas 
the outer is the systematic error added in quadrature. The relative uncertainties from the 
fragmentation model are added in quadrature to the systematic uncertainty of the data, 
whereas all other extrapolation uncertainties are indicated by a shaded band at the bottom 
of a certain Q 2 bin of F% as shown in figure [3^) and b). The prediction of HVQDIS using 
the proton PDF CTEQ5f3 is shown in Fig. [3^,) by the dark solid line and describes the data 
rather well. In addition also MRST04FF3 was used with HVQDIS, which underestimates 
the data slightly at low x and Q 2 . In case of CASCADE the unintegrated gluon distribution 
AO was used and is indicated in Fig. [5]d) by the solid line, which describes the data nicely. 
A comparison of F£ using HVQDIS or CASCADE shows differences dominantly at high x 
and especially a steeper slope towards low x for the values extracted with HVQDIS. 

4 Conclusions 

The charm contribution, i 7 ^ , to the proton structure was determined using the measured 
cross sections of D* production in y and Q 2 . Two theoretical models, HVQDIS and CAS- 
CADE, were used for the extrapolation in pr(D*) and r/(D*) of the visible cross section 
to the full phase space. Differences between the extrapolation factors of these models are 
observed at high x. The uncertainty of the fragmentation model has been estimated from 
results of an HI measurement. The total theory uncertainties are much larger than the ex- 
perimental errors. Although especially at low x the predictions from HVQDIS with different 
proton PDFs differ slightly the discriminating power of the current analysis is not sufficient. 
Otherwise the comparison of the data to the predictions of HVQDIS and CASCADE shows 
a good agreement. 
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Figure 3: F% as determined from the inclusive D*^ meson cross sections using HVQDIS a) or 
CASCADE b). The error bars refer to the statistical (inner) and systematical error (outer) 
added in quadrature. The solid lines represent the predictions of F% using appropriate 
parton densities for HVQDIS and CASCADE. The shaded band shows the relative change 
of the central value due to variation of the model parameters. 
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